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TECHNICAL NOTE 339O 


SECOND-CKDER SUBSONIC AIRFOIL -SECTION THEORY 
AND ITS PRACTICAL APPLICATION 
By Milton D. Van Dyke 


SUMMARY 


Several recent advances in plane suLsonic flow theory are combined 
into a unified second-order theory for airfoils of arbitrary shape. The 
solution is reached in three steps: The incompressible result is found 

by integration, it is converted into the corresponding subsonic com- 
pressible result by means of the second-order compressibility rule, and 
it is rendered uniformly valid near stagnation points by further simple 
rules. Solutions for a number of airfoils are given and are compared 
with the results of other theories and of experiment. A straightforward 
conputing scheme is outlined for calculating the pressures on any airfoil 
at any angle of attack. 


INTRODUCTION 


Thin-alrfoil theory provides a useful first approximation to the 
incompressible flow past two-dimensional airfoils, and the results can 
be immediately extended to subsonic compressible flow by the Prandtl- 
Glauert rule. It is natural to attempt to improve this simple theory by 
successive approximations so as to increase its accuracy for thicker air- 
foils and higher subsonic Mach numbers. There results a series expansion 
of the flow quantities in powers (supplemented in some cases by logarithms) 
of the airfoil thickness ratio, camber ratio, and angle of attack. 

For incompressible flow, the higher-order theory has been studied by 
various writers, in particular Riegels and Wlttlch (refs. 1 and 2 ) and 
Keune (ref. 3) . A less straightforward series of approximations was 
developed by Goldstein (ref. ^<-) . Perhaps the most concise exposition of 
higher-order incompressible thin-alrfoil theory is given by Lighthlll 
(ref. 5). 

For subsonic compressible flow, the corresponding analysis was first 
undertaken by Gortler (ref. 6), followed by Hantzsche and Wendt (refs. J 
and 8), Schmieden and Kawalki (ref. 9), Kaplan (refs. 10 and ll) , and 
Imal and Oyama (refs. 12 and I3) These investigators treated only spe- 
cffic Bimple shapes by rather laborious analysis. Later, it was discovered 

1 These historical references are intended to be representative rather 
than eadiaustive. 
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that particular Integrals of the second-order Iteration equation can he 
expressed in terms of the first approximation (refs. l4 and 15)* This 
permits the second-order suhsonic solution for any profile to he given 
in terms of integrals (refs. 15 and l6) . However, the resulting solutions 
are incorrect everywhere for airfoils with sta^iation points, for reasons 
to he discussed later. 

Recently Hayes (ref. 17 ), improving on a result of Imal (ref. l8), 
has given a second-order similarity rule for surface pressure that implies 
a second-order extension of the Prandtl-Glauert rule (ref. 19) • This 
remarkable result was overlooked hy earlier investigators because they 
did not calculate surface pressures, hut were content with finding surface 
speeds, for which the second-order compressibility rule is more compli- 
cated. These rules reduce the second-order problem of suhsonic compress- 
ible flow past airfoils to the corresponding incompressible problem. 

However, the solution by successive approximations breaks down near 
leading and trailing edges if there are stagnation points. The result is 
therefore merely a formal series expansion, which fails to converge near 
the edges. Ifi first-order theory spurious singularities arise at stagna- 
tion edges, but it is known how they can be taken into account, since 
they are integrable. In the second approximation, however, these singu- 
larities are Intensified, so that at round edges they are no longer inte- 
grable. In any case, the calculated speeds and pressures are incorrect 
near such edges > and more so in the second approximation than the first. 
Moreover, in subsonic compressible flow the second approximation may be 
Incorrect everywhere as a consequence of the defects in the first approx- 
imation. 

For round edges in Incompressible flow, previous Investigators have 
shown how these defects can be corrected. Rlegels (ref. 2) gave a simple 
rule that renders the first-order thln-airfoll solution valid near the 
edge. Llghthill (ref. 5) gave an equivalent rule for the second approx- 
imation. Recently, corresponding rules have been developed for higher 
approximations, for sharp as well as round edges, and for subsonic com- 
pressible flow (ref. 20) . 

It is the aim of this paper to combine these recent advances into a 
unified theory. There results a uniform second approximation to subsonic 
flow past any profile at angle of attack, expressed in terms of integrals 
that can, if necessary, be evaluated numerically. It may be noted that 
the resulting solution is now generally believed to be valid only below 
the critical Mach number - that is, for purely subsonic flows. Although 
only flow quantities at the airfoil surface are considered here in detail, 
the entire flow field can be treated in the same way. 

For numerical ccanputation, the method initiated by Germain (ref. 21), 
and extended by Watson (ref. 22), Thwaites (ref. 23), and Weber (ref. 24) 
appears to be the most useful. It requires a knowledge only of the air- 
foil ordinates at a specified set of points. A straightforward scheme, 
based on this method, is given for computing the second-order subsonic 
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solution for any airfoil. The reader interested only in calculating a 
specific case, without necessarily \xnderstanding the theory, can turn 
directly to the section "PRACTICAL NUMERICAL COMPUTATION" on page I9. 


THEORY 


From the preceding remarks it is clear that the solution is reached 
in three steps. First, the formal second-order Incompresslhle solution 
is found hy integration. Second, this is converted into the correspond- 
ing subsonic compresslhle solution hy means of the second-order ccanpress- 
ihility rule. Third, this is modified near stagnation points hy the 
appropriate rules for round or sharp edges. These three steps will he 
considered successively. 


Formal Incompressihle Solution 


The expansion of the velocity components in a formal series of powers 
of the airfoil thickness ratio, camher ratio, and angle of attack has heen 
discussed in detail hy Lighthill (ref. 5) . It will suffice here to sum- 
marize his results for the second approximation. We mainly follow his 
notation except to make it more mnemonic, and to suppress- his parameter 
e characteristic of the airfoil thickness, which is only convenient in 
the detailed analysis. 


Accordingly, consider an airfoil of moderate thickness and camher at 
a moderate angle of attack to a uniform subsonic stream (sketch (a)). It 
is essential that the x axis he 
chosen to pass throu^ both the lead- 
ing and trailing edges. Let the upper 
and lower svir faces of the airfoil he 
described hy 

y = Y(x) = C(x) ± T(x) ( 1 ) 



where c(x) describes the mean camher line and T(x) the thickness. The 
airfoil extends over the Interval A<x<B, which is usually conveniently 
taken to he either -l<x<l or 0 <x<l. All symbols are defined in 
Appendix A. 


First-order solution .- In the first approximation of thin-airfoll 
theory, the condition of tangent flow at the alrfpll surface is imposed 
on the two sides of the chord line y = 0 rather than at the smrface, and 
requires that 


Zl 

U 


y=o 


Y'(x) = C'(x) ± T'(x) 


( 2 ) 
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The corresponding horizontal velocity disturhance on the chord line, which 
is required for calculating the surface pressure, consists of a term 
associated with the airfoil thickness, and another associated with its 
camher and angle of attack. For the thickness 


uit _ 1 T*(|) dg 

U - « X - I 

and for the camher and angle of attack 

^ - xV ^^ r 1 r® / 'l - aV ^ ^ C»(|) d| 


( 3 ) 

(^) 


The latter result la due to Munk (ref. 25 ) and the former was apparently 
first given hy Squire in a paper that is still not generally available. 
Cauchy principal values are indicated in each integral. 


The surface speed Is then given to a first approximation hy 

= 1 + ^ 
u U “ u 


( 5 ) 


Second -order solution .- In the second approximation, the tangency 
condition is transferred from the airfoil surface to the. chord line hy 
Taylor series expansion. The condition on the second-order increment In 
vertical velocity is thus found to he 



C’aCx) ± T’sCx) 


(6a) 


where 


Ca(x) 


^it 

U 


C + 


uic 

U 


T 




> 


Ta(x) = ^ T 


Ulc 

U 


(6h) 


(We depart here from Llghthill’s notation in order to emphasize that the 
functions Cg and Tg are effectively the camher and thickness for some 
fictitious airfoil.) The problem is identical with that in first-order 
theory except for the condition at infinity, which is readily disposed 
of. Thus, corresponding to Tg is the increment in horizontal velocity 

^ - i Tp(l) d| 1 

U « X - I ■ 2 


(7) 
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and corresponding to 

_ 1 'Z® A - AV ^g C*g(g) dg 

U - It Vx - Ay' \ ^ ^ ^ ^ 

The velocity components on the surface of the airfoil include also 
terms arising from the transfer from the chord line to the surface, which 
is again effected hy Taylor series expansion. Hence the surface speed is 
given to a second approximation hy 


£2 = 
U 


l+^±^ + ^±^+(C±T) (C ± T") 


\ (c ± 


T')2 (g) 


Airfoil integrals .- The Inccmpresslhle solution to second order (or, 
indeed, to any order) is thus reduced to a succession of "airfoil inte- 
grals" typified hy equations (3)^ (?)» and (8). Goldstein (ref. 26) 

emphasizes that in first-order theory these Integrals can he evaluated 
analytically for practically every profile for which formulas have ever 
heen proposed. In second-order theory this appears to he true to a some- 
what lesser extent, although the labor of calculation becomes great except 
for simple shapes. Often the integrals are most readily evaluated hy 
guessing (u - iv) as a function of the complex variable (x + iy) that has 
the required behavior on the chord line. A short table of airfoil inte- 
grals useful for finding second-order solutions is given in Appendix B. 
Others can he found in references 26 and 27- 


For complicated profiles, exact analytic evaluation of the integrals 
may he impossible or excessively laborious. Then mmierlcal integration 
may he resorted to, or the profile can he approximated hy a simpler shape 
that can he treated analytically. The most useful numerical procedure is 
apparently that originated hy Germain and simplified and extended hy 
Watson, Thwaltes, and Weber, In this method the airfoil ordinates are 
approximated hy the trigonometric polynomial 


N-i 


Y “ Co + I (cpcos r© + trsin r©) + 


r=i 


cj^cos W© 


where © is the angle indicated in 
sketch (h) . The coefficients cr (f 
camber) and tr (for thickness) are 
chosen to give the actual ordinates at the 
2N points for which © = mit/N. In this 
way it is found that the airfoil integrals 
can he expressed approximately as sums of 
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the airfoil ordinates at 'certain pivotal points multiplied by standard 
inf luence -coefficients . The details of this method, as adapted to thin- 
airfoll theory, are given in Appendix G. The numerical computing pro- 
cedure is outlined in the last section of this paper. 


Second-Order Compressibility Rule 


The second-order counterpart of the PraSldtl-Glauert compressibility 
rule is implicit, in an extension of transonic similitude that was initi- 
ated by Imal (ref. 18) and carried to completion by Hayes (ref. 17 ) . 

Imal sought to improve the transonic similarity rule by retaining in its 
derivation all terms proportional to the square of the airfoil thickness 
except one appearing in the condition of tangent flow at the surface. 

The correlation of experimental data was not appreciably improved, which 
led him to suggest that the neglected second-power term should also be 
included. This probably cannot be done. However, in attempting merely 
to reproduce Imai's result as announced before pu'blication, Hayes actually 
included that term in a second-order rule fnr surface pressure. 

Hayes’ result is that for two-dimensional subsonic or supersonic flow 
the ratio of the second-order to first-order pressure teim on the surface 
is proportional to the parameter 


T 

(|l - 


— g ^ + 2(1 - m 2 ) 


(11) 


where r is some measure of the thickness, camber, or angle of attack. 
Now at subsonic speeds the first-order pressure term is related to its 
value in incompressible flow by the Prandtl-Glauert rule. Combining 
these two results yields the second-order compressibility rule (ref. 19) . 

In incompressible flow the second-order surface -pres sure coefficient 
has the form 


Cpo(x) = Cpi(x) + ACpg(x) (l 2 a) 

where the first-order term contains linear terms in thickness, cam- 

ber, and angle of attack, and the second-order increment ACpg contains 
their squares and products. Then for the same airfoil in subsonic com- 
pressible flow, according to the compressibility rule, the pressure coef- 
ficient is 


Cpjj = KiCpj^ + K2(ACp2) 


(12b) 


where 
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Ki = 


y/l - m2 P 


K= 


_ (y + 1)M^+ 


( 12 c) 


It has been pointed out that the formal thin-airfoil series requires 
modification near stagnation edges. The modification must in general he 
performed on the speed rather than the pressure. Hence the compressibil- 
ity rule for surface speed is required. It is readily found frcxn the 
above rule for pressorre by considering the small-disturbance series form 
of Bernoulli's equation for compressible flow. Thus it is found that if 
the surface speed ratio in incompressible flow is 


= 1 + ^ ^ 

U U U 


then at subsonic speeds 




with 

Kg - 1 _ m2 (7 + 1)M^ + 43^ 
2 83* 


(13a) 


(ist) 


(13c) 


This rule is seen to lack the fundamental simplicity of the rule for pres- 
sure . 


Modification for Stagnation Edges 


Thin-airfoil theory is known to fail near leading and trailing edges 
if there is a stagnation point. The flow is actually brought to rest, but 
thin-airfoil theory predicts infinite speeds instead. If r is the dis- 
tance from the edge, the velocity contains powers of for a round 

edge and for any leading edge with. flow arcund it (associated with angle 
of attack), and powers of In r for a sharp edge. First-order theory 
contains first powers of these singularities, second-order theory their 
squares, and so on, so that the formal thin-airfoil series diverges in 
some neighborhood of the edge. Hot only are the velocities and pressure 
incorrect near stagnation edges, but nonlntegrable singularities appear 
in the higher-order expressions for aerodynamic forces. 
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False sulasonlc. eolutlons .- Even more serious difficulties may arise 
in subsonic compressible flow, where the infection spreads in some cases 
BO that the formal second-order solution is Incorrect not only near the 
edges but over the entire airfoil surface. Thus, using the particular 
integral of reference l 4 . Harder and Klunker gave an expression for the 
second-order solution for any symmetric airfoil at zero angle of attack 
(ref. 16). However, they noted that their expression does not apply to 
round-edged airfoils, for which it contains divergent integrals. A more 
deceptive defect appears if their. express ioia is applied to a sharp-edged 
airfoil such as a biconvex section; then the predicted surface speeds are 
finite (except near the edges) but incorrect everywhere by a term propor- 
tional to M^. This defect arises from the fact that near the edges the 
first-order source_distribution is not approximately the airfoil slope, 
as is assumed in thin-airfoil theory. The second-order solution involves 
the derivative of the so\n:ce strength which, as indicated in sketch ( c) , 

has sharp peaks that are missed by 
thin-airfoil theory. It is enough to 
take account of this shortcaning in 
even the crudest fashion. Thus, if 
the region of integration is extended 
an infinitesimal distance beyond the 
edges to include the pulses (Dirac 
delta functions) of the thin-airfoil 
approximation. Harder and Klunker 's 
expression yields a solution that is 
correct to second order except in the 
vicinity of the edges. 

Keune has discovered an alterna- 
tive particular integral containing 
the stream function rather than the 
velocity potential, and so has obtained another expression for the second- 
order solution (ref. I5) . Because the tangency condition is one degree 
smoother for the stream function than the 'yelocity potential, his expres- 
sion yields the correct result (except near stagnation edges) for sharp- 
edged shapes. It fails, however, for round-edged shapes, so that his 
solution for subsonic flow past an ellipse is incorrect everywhere. 

Both these expressions can be manipulated by partial integration so 
as to be correct except near stagnation edges. However, the result is 
simply that obtained by applying the second-order ccanpressibility rule to 
the expressions for second-order incompressible flow. Hence these more 
serious difficulties are of no further concern here. They do serve, how- 
ever, to warn of the danger of false second-order solutions in more com- 
plicated problems . 

Modification for incompressible flow .- For round edges in incompress- 
ible flow, Riegels (ref. 2) and Lighthili (ref. 5) have given simple rules 
that render the formal thin-airfoil solution uniformly valid. The result 


Actual values 


Source 

strength 


00 


00 


Its 

derivative 


Thin -air fait 
approximation 


Sketch (c) 
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is at least a first approximation to the flow distiirbances near the edge, 
Riegels found a rule for first-order theory hy considering the conformal 
mapping, and Llghthlll found a rule for second-order theory hy considering 
a contraction of abscissas that shifts the thln-airfoil solution by half 
the radius of the edge. 


In reference 20 corresponding rules have been developed for higher 
approximations, sharp edges, and subsonic compressible flows. The tech- 
nique used there is to consider the exact solution for some simple shape 
that approximates the airfoil in the vicinity of its edge. The ratio of 
the exact solution for the simple shape to its formal thin-airfoll series 
expansion serves as a multiplicative factor that corrects the series 
expansion for the actual airfoil. The result should then be simplified 
Insofar as possible. The relevant rules will be summarized here; the 
details are given in reference 20. 


A round-nosed airfoil can be closely approximated by a parabola whose 
axis coincides with the initial camber line (sketch (d)). The exact solu- 
tion for Incompressible flow past the 
parabola (resolved into streaming and 
circulatory components) leads to the 
following rule that converts the formal 
second-order solution "q^ for siirface 
speed into a uniformly valid approxima- 
tion 


U 


(; 


Xn ± \ </2pXn 


Xq + p/2 ± A T^Xo, 




Here Xq is the abscissa measured from 
the edge into the airfoil, p is the edge 
radius. Tv is the initial angle of cam- 
ber, and the ±.. signs refer, as usual, 
to the upper and lower surfaces. 



Sketch (d) 


This mile yields a uniform second approximation to the disturbances 
everywhere (except at the other edge, where edditional modification may 
be required) if the rate of change of curvature of the profile is con- 
tinuous, which means that near the edge the thickness has the form 


T(x) = bi + baXo®'' ^ + . . . (15) 

However, airfoils of the NACA four- and five-digit series violate this 
requirement, their leading edges having the Initial form 


T(x) = bi + bsXo + t>3Xo3/2 + 


( 16 ) 
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Hence the rule yields only a first approximation near the edge (while , 
leaving a second approximation elsewhere) and can therefore be replaced 
by the simpler form 


£2 _ 
U 


xo 


Xq + p/2 


,1/ 2 




(17) 


which is Lighthill' a rule. 

The rule for first-order theory is obtained by dropping the term 
p/4xq from equation (iT) • However, it is then advantageous to use the 
alternative form due to Riegels (ref. 2), which is correct to the same, 
order but a great deal more accurate (see ref. 20) . It is simply 


^ = (cos ^) ^ 


(18) 


where t] is the angle of the airfoil surface. Analogous alternative 
rules can be found for the second-order theory, but their shortcomings 
probably outweigh their slight advantages of accuracy and simplicity, so 
they will not be considered here . 


The modification for a sharp edge is found by considering incompress- 
ible flow in an angle. If the edge is a trailing edge with Kutta condi- 
tion enforced, or a leading edge at the ideal angle of attack, so that 
there is no flow around it, the second -order rule is 


£2 = 

U 


6 


Xo 




-£2 _£ _£±. In 
U n u 


Xo 



Xo 


1 

2 


In^o 


(19) 


where 8 is the semi vertex angle, and "qj! is the first-order solution. 
Otherwise, the circulatory part of the formal second-order solution, which 
consists of the terms singular like must first be corrected sep- 

arately by the rule 


£2 = 
U 


Xo 


5 

2(jt-5)' 


R 11^11 

_£2 _ A 

U 2it U 


In Xo 


( 20 ) 


after which the remainder is corrected by equation (l9)* 

Airfoils with two stagnation edges can be treated either by applying 
the appropriate correction separately at each edge, or by combining the 
rules. The ccanbined rule for two round edges is given in equation (24) 
of reference 20. Similarly, for a round edge at x = -1 and a sharp edge 
(with Kutta condition) at x = 1, the combined rule is 
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(1 - xj-'-S 


1 4- X ± «/gp(l + x) 

1 + X + p/2 ± A j/2pCl + x) 


r ft 

J _a2.° _2i 
\ u It V 


p 

4(1 + x) 


— i ln^(l - x) - In (1 - x) |- 


For two sharp edges of equal angle, hoth with Kutta condition (as for a 
biconvex airfoil at zero angle of attack), and located at x = ±1, the 
combined rule has the form of equation ( 19 ) with Xq replaced by 


If the pressure is required, it must be calculated from the speed 
using the full Bernoulli equation because the disturbances are no longer 
assumed to be small. 

Modification for subsonic flow .- For round noses, these rules are 
extended to subsonic speeds by considering ccmpressible flow past a parab- 
ola. Thus the counterpart of equation (l4) is found to be^ 


U ^ V p ^ J L u 


U 


i Ko 


Kx^ - 


Here, as indicated in sketch (e) , 

Q is the speed ratio on the sur- 
face of a parabola in a \iniform 
subsonic stream of Mach nijmber M, 
with circulatory flow proportional 
to a. The combination 
Xq = Xo ± t/^p^o appears in the 
first argument of Q, as it does 
in equation (l4), because of the 
connection indicated in sketch (e) 
between abscissas of a surface 
point measured along the x axis 
and along the axis of the parabola. 
(See also eq. ( 9 ) of ref. 20.) The 
factor a is proportional to the 


/ 





Sketch (e) 


^It should be noted that as M -» 0 this rule reduces not to the 
incompressible rule of equation (l4), but to an alternative that is 
entirely equivalent up to terms of second order. See footnote 7 of ref- 
erence 20 . 
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angle of attack of the actual airfoil measured from the ideal angle at 
which the stagnation point coincides with the vertex. It must he found 
as the coefficient of in the first-order solution ”qi/u. 

The function Q is not known exactly, hut a satisfactory approxima- 
tion is given hy the Janzen-Rayleigh solution in powers of M®. For the 

special case of corresponding to the ideal angle of attack 

where a = 0 , the solution has heen calculated to order M* hy Imai 
(ref. 28 ) . It is tabulated briefly in reference 20 for y = where 

it is denoted hy Q(x/p, M) . For other angles of attack, the function 
Q to order can he extracted hy a limiting process from Kaplan's 

solution for an inclined ellipse (ref. 29 ), which gives, with Xq/p » cp, 

= T)f 


0((p,^k,M) = -=^== 

Vl + 2sp 


^ " 2(1“+ ' ^T ^ 


1 + 

i + li<p2 


^ - I Ilr + cpt) In - (W - 1 


2t«/2r^ tan“^ »/2t J- 


( 23 ) 


The rule for sharp noses in subsonic flow can he found hy considering 
compresslhle flow in an angle. However, this basic solution is not yet 
available. For practical purposes the correction is probably negligible 
since it is appreciable over a much smaller neighborhood of a sharp edge 
than a round one. Moreover, sharp edges are usually trailing edges, in 
which case the details of the flow are altered by viscous effects. 


EXAMPLES: CCMPARISOH WITH EXPERIMENT 

AND OTHER THEORIES 

Inccaapresslble Flow 


It has been seen that the solution for subsonic flow depends on that 
for incompressible flow. It is therefore pertinent to test the second- 
order theory in the case of incompressible flow, where it can be checked 
against the exact results of conformal mapping. 
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Ellipse . - Consider an ellipse of thickness ratio t with the Inter- 
val -1 < X < 1 as chord line. It is described hy 


= ± T - x2, -1 ^ X < 1 


(2lf) 


Suppose that the Kutta condition la satisfied - the rear stagnation point 
coincides with the end of the major axis. Then the first-order solution 
for surface speed is found, from equations ( 3 ), (4), and ( 5 ), together 
with Appendix B, to he 




- X 


+ X 


(25a) 


Proceeding with equations (6) to ( 9 ) gives the formal second-order result 


^ = 1 + T ± 


X 1 X2 

— T2 ± OtT 


+ X 2 1 - X* 


vfil - 1 (25b) 


This can he checked hy expanding the exact result, which is 


= Cl + t) < \/1 - x2 cos g ± (1 - x) sin a 

U ^1 - x2 + t2x2 


(26) 


The formal second-order solution clearly breaks down near the ends of the 
ellipse. It is converted into a xiniformly valid second approximation hy 
applying equation (l4) twice in succession, or using the combined rule of 
equation (24) of reference 20 , which gives 


U 



1 - x2 
- x2 + 


t2 


(1 + 


t) 



These approximations are conpared in figure 1 with the exact solution 
for an 18 -percent-thick ellipse (which has nearly the same nose radius as 
an NACA 0012 airfoil) at zero angle of attack. The precipitate descent 
of the formal second-order solution toward negative infinity is just dis- 
cernible near the nose and is eliminated in the modified theory. It 
should he mentioned that the first-order theory modified according to 
Rlegels' rule (eq. ( 18 )) happens to give the exact result for an ellipse. 



14 


WACA TN 3390 


Symmetrical Joiikovskl airfoil .- To seccmd as veil as first order a 
symmetrical Joiikovskl airfoil of thickness ratio t 1 b described by 


y = + Ti(i - x) x2, - 1 < X < 1 

> 

Ti = JL T = O.7698T 


(27) 


By the foregoing procedure, the formal second-order solution is found to 
be 


¥ " 1 + T^l - ac) ± 


1 

2 


Tl2 


1 - X 
1 + X 


(1 + 2x)^ + 


Srxctxy^-^ - I 0,2 ( 28 a) 


where the first three terms give the first-order solution. Modifying this 
according to equation (l4) with Xq = 1 + x and p = 4 tx 2 (and A = O) 
gives the uniformly valid second approximation 


I2 - / ^ ^ 

U yi + X + &T 3.2 


1 + 


Ti(l - 2x) + I 




2x)2 + 





(28b) 


In figure 2 these approximations are compared with the exact solution 
(ref. 30 ) for a 12-percent-thick section at zero angle of attack. The 
effect of the modification on the second-order result is not discernible 
to this scale. 


Biconvex airfoil .- To second order a symmetrical biconvex airfoil of 
thickness ratio t bounded by either circular or parabolic arcs is 
described by 


y s= ± t (1 - x 2 ), -1 < X < 1 


(29) 


The formal second-order solution is found to be 
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^ = X + I r(2 - X (l 


^(s.xln^ 


1 - m2 1±J£ 

jr2 1 - X 


- (1 - x2) 



(l + 2 x) ln ^'^^-4 

J 1. + X 

1 - X 

^ 0" 


. ia2 
2 


(30a) 


In deducing from this a miformly valid second approximation, the terms 
independent of a are treated by the rule for combined equally sharp 
edges that was described just after equation ( 2 l), with 6 = 2t. The 
terms in a are modified according to equation (20) with xq = 1 + x. 
(Notice that no modification of these terms is required at the trailing 
edge.) The result is 


2T 


U 


^ = (1 - x2)” + [2 - (i + x) m (i + x) - (1 -x) m (1 -x)] + 




3 - (1 - x2) - 3(1 + x) In (1 + x) - 3(1 - x) m (1 - x) + 


r (1 + 3x)(i + x) m2 (1 + x) + i (1 - 3x)(i - x) m2 (1 - x) + 
4 4 


I (1 - x2) m (1 + x) m (1 - x) 


I 

^[2(i + x) m(i + x) -(1+ 2x) m (1 - x) - 4 ] 1 - 


(1 + x) 


«-2T |l - 


(30b) 


These approximations are compared m figure 3 with the exact solution 
(ref. 31) for a circular-arc airfoil 18 percent thick at zero angle of 
attack. Although the vertex angles are large in this example, the modifi- 
cation of the second-order solution is appreciable in such a small neigh- 
borhood of the edge that it would be invisible even on a much larger plot. 

NACA OOXX airfoils.- Symmetrical airfoils of the NACA OQXX family 
(such as the NACA 0012 ) are nationally defined for the Interval 0 < x < 1 . 
The airfoil of thickness ratio r is described by (ref. 32) 


y = ± T(x) = ± r(bi + bgx + b4x2 + bgx^ + bax*) , 


0 < X < 1 (31) 
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where 

hi = l.h8A50 

bg = -0.63000 
b4 = -1.25800 
bg = 1.42150 
be = -0.50750 

With the aid of Appendix B the first-order solution is foimd to be 


u ^ It 


- T’(x) In 

T 



1 + <y/x 

7 x 




(3bg + 2ba)x - 4b^2^ ± 0. ^^ - " - — (32a) 

in agreement with the result given by Goldstein (ref. 26 ) . Applying 
Riegels’ rule (eg. (18)) renders this a uniformly valid first approxima- 
tion except very near the trailing edge. 

The second-order terms in thickness, in addition to being very com- 
plicated, involve integrals that apparently cannot be evaluated in terms 
of tabulated functions. Accordingly, the second-order terms have been 
calculated using the Germain -Watson-Thwaites -Weber numerical method dis- 
cussed in Appendix C, with N = I6. The accuracy of this approximation 
is assured by the fact that cruder approximations modify the numerical 
results only slightly, as will be seen in a later example. 

The formal second-order solution for surface speed therefore has the 

form, 

^ = 1 + tQt ± a + r^r ± | ( 32 h) 

where values of Qt from equation (32a) and approximate numerical values 
of Qtt Q^a 


X 

Qt 

Qtt 

Qra 

X 

Qt 

Qtt 

Qra 

0.025 

1.943 

-9.00 

8.80 

0.50 

0.900 

-0.135 

0.32 

.05 

1.836 

-3.35 

5.55 

.60 

.697 

-.220 

.11 

.10 

1.714 

- 1.00 

3.25 

.70 

.485 

-.315 

-.08 

.20 

1.510 

-.090 

1.65 

.80 

.238 

-.410 

-.23 

.30 

1.309 

.010 

1.00 

.90 

-.124 

-.420 

-.33 

.40 

1.106 

-.060 

.58 

.95 

-.440 

-.360 

-.34 
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Applying equation (l4) with p = I.IOI 87 8 « 1.16925 t (and 
A = 0 ) yields a unifoimly valid approximation. However, as discussed 
previously, the curvatvire of the profile does not vary continuously near 
its nose, so the result is only a first approximation there, though a 
second approximation elsewhere. 

The various approximations are compared in figure 4 with the result 
of a "long and elaborate calculation" by conformal mapping for the 
NACA 0012 airfoil that is given by Goldstein (ref. 33 ) . Again the effect 
of modifying the second-order solution is indiscernible. Also shown is 
the "exact" solution tabulated in reference 3^» The agreement between the 
first-order solution with Rlegels' rule, the second-order solutions, and 
Goldstein’s calculation leaves little doubt that his is the more accurate 
of the two "exact" solutions. 


Compressible Plow 


When extended to subsonic compressible flow, the preceding examples 
can all be compared with other theories or with experiment. As before, 
the comparisons will, for simplicity, be made only for zero angle of 
attack. 

Ellipse . - Applying the second-order compressibility rule of equa- 
tions (I 3 ) to the Incompressible solution of equation (25b) gives as the 
formal second-order solution for the speed on an elliptic cylinder 


_i2 

U 


1 + Ki( 





+ 


ra( 2K2 


i) 


yrr 


- X 

X 


- a,‘ 


1 + (2K:, - l)x 
2(1 + x) 


( 33 ) 


For zero angle of attack the maximimi speedy occurring at midchord, is 
given by 


^ A 


= 1 + K]_t + 


max 


i (Ks - 1)t2 


(34) 


in agreement with the result of Hantzsche and Wendt (ref. 7 ) • Hantzsche 
has also calculated the third-order solution for the maximum speed at zero 
angle of attack (ref. 8) . Values of the maximum speed ratio calculated 
from these and other approximations for a 10-percent-thick ellipse at zero 
angle of attack are 
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M = 0.70 

M = 0.75 

M = 0.80 

First-order theory (or Prandtl- 




Glauert rule applied to exact 




incompressible value of l.lOO) 

l.l40 

1.151 

1.167 

Karm^-Tsien rule 

l.ll)5 

1.166 

1.184 

Second-order theory 

1.151 

1.166 

1.189 

Third -order theory 

1.151 

1.172 

1.198 


Here the Karman-Tslen rule has been applied to the exact incompressible 
value of the pressure coefficient, and the speed ratio then calculated 
from Bernoulli’s equation. It is to be anticipated that second-order 
theory is more accurate than any of the compressibility correction for- 
mulas such as the Karman-Tsien rule, because it allows for a dependence 
on the particular airfoil shape and on the value of 7. This is seen to 
be true for the ellipse. 

In the same way the second-order solutions are readily calculated 
for the Joukowskl and the biconvex airfoils, and are found to agree with 
the results that Hantzsche and Wendt obtained by laborious analysis. 

NACA 0012 airfoil .- The formal first- and second-order solutions for 
NACA OOXX airfoils in subsonic flow are easily obtained from equations (I3) 
and (32) . The second-order solution can then be rendered uniformly valid 
near the nose using equation ( 22) , although again the modification is sig- 
nificant in only a very small regicm of the nose. 

For the KACA 0012 airfoil at zero angle of attack, Emmons has calcu- 
lated the flow field at Mach numbers of 0, O.TO, and 0.75 using the numer- 
ical relaxation method (ref. 35) • last of these Mach numbers is super- 

critical, so that the flow contains shock waves, and is beyond the scope 
of the present theory. The pressure distribution calculated by the relax- 
ation method for M = 0.70 is compared in figure 5 with the results of 
first- and second-order theory and various other approximations. The 
relaxation solution for incompressible flow is also shown in comparison 
with Goldstein' s "exact" solution, and is seen to be inaccurate near the 
nose. The solution for M = O.7O probably contains similar inaccuracies, 
however, just as for the ellipse the pressure coefficients calculated by 
second-order theory may be slightly less negative than the true values 
near their minimum. 

Experiments on NACA OOI5 airfoil .- Experimental pressure distribu- 
tions in two-dimensional flow over the NACA 0015 airfoil at high subsonic 
speeds are reported in reference 36. For zero angle of attack, the criti- 
cal Mach number is approximately O.7O. The measurements at this Mach num- 
ber are compared in figure 6 with the results of first- and second-order 
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theory and of the two common compressihility correction formulas applied 
to the Incompressible flow values tabulated in reference 3 ^* Unfortu- 
nately, the model was imperfectly constructed, and the ordinates were 
inaccurate near the nose and midchord. Otherwise, the measured pressiares 
are in satisfactory accord with either second-order theory or the results 
of the Karman-Tsien rule. 


HIACTICAL NUMERICAL COMPUTATION 


The following computing procedure yields the formal second-order 
subsonic solution for the surface speed or pressure on any airfoil at any 
angle of attack. It requires a knowledge only of the airfoil ordinates 
at seven points along the chord. It is based on the foregoing theory 
together with the ntunerical method of Germain, Watson, Thwaites, and Weber 
that is discussed in Appendix C. 


Computing Procedure (N = 8 ) 


(1) Tabulate the ordinates Yu and Yi of the upper and lower sur- 
faces at the seven pivotal points xn listed in table I. (The x axis 
must pass through the leading and trailing edges.) 

(2) Calculate the corresponding values of 


T = I (Yu - Yz), C = I (Yu + Yz) 


(35) 


(3) Using the influence coefficients of tables II and III, calculate 


¥= I 




s=i 


-I 


S=1 


■£ 




dnsCf 


ensTs C* = > fnsCs 

B=1 S=1 




(36) 


r T 

T" = ^ Shs'I's C" = huBCs 


s=i 


s=i 
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4. Using table I (with a in radians), calculate 


^ic _ ^ic* 
U ~ U 


+ a 


- X 


(37) 


5 . Calculate 


Tj 


^ T + ^ C, 
U U 


Co = c + T 
^ u u 


(38) 


6. Using the influence coefficients of _ tables II and III, calculate 


Ugt 

U 


f T 

= ^ CnsT'ss ■ ^ ^ns*^ 

S=1 S=1 


ss 


7. Using the compressibility factors of table IV, calculate 
•^ = 1 + Ki ^ + CC + TT" + ^C'^ + + 

+ CT" + 




urt 

U 




Kl ^ + Ks ^ 


(39) 


C"T + C’T' ) + (Ka - 1) 


^it ^ic 
U U 


(ho) 


or 


’°P2 - 


2 ^ + 2 ^ + 2(C ± T)(C" ±-T") 


(C ± T')' 


Alt 





(hi) 


The ± signs refer to the upper and lower surfaces of the airfoil. 
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Remarks 


1 . The summations of steps (3) and (6) are conveniently carried out 
hy tabulating C and T in columns that can he matched with successive 
columns of tables II and III while cumulative multiplication is carried 
out on a desk calculating machine. 

2 . If the airfoil slope and second derivative can be found more 
directly, step (3) can be simplified by omitting the calculation of T', 
C, T", and C". ‘ 

3. The results near the nose can be rendered valid by the use of 
the modification of equation ( 22 ), but the effect is often insignificant. 

ij-. Seven pivotal points yield sufficient accuracy for most purposes. 
If conditions near the nose of a thin airfoil are of interest it may be 
necessary to repeat the computation using I5 pivotal points. (The values 
of T and C already calculated can be used again.) Table I gives the 
additional pivotal points and angle -of -attack solution, and table V gives 
the Influence coefficients for calculating u/Uj additional tables can be 
prepared for calculating the airfoil slope and curvatvire if they are not 
known otherwise. 

5. The above scheme is designed for calculating a single case. If 
the same airfoil is to be calculated at more than two argles of attack, 
it is economical to subdivide the computing scheme to separate terms in 
a and oc?. Similarly, the scheme should be subdivided if more than two 
thickness or camber ratios are to be calculated for the same family of 
airfoils . 

6 . For MCA airfoils T is the basic thickness and C the camber 
line. To second order it is Immaterial that the thickness is added normal 
to the camber line rather than to the chord line. 


Example 


The following table gives the complete computing sheet for calculat- 
ing the first- and second-order increments in surface speed for an 
MCA OOXX airfoil (of unit thickness ratio) at zero angle of attack and 
zero Mach number: 
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n 

X 

T 

^it _ Aqi 

TT." “IT 

T* 


Ta 

uat 

nr 

Aqa 

nr 

1 

'S 

00 

CT) 

0 

• 

0 

0.26316 

1.8936 

3-0914 

-54.721 

0.49832 

4 . 8 o 4 o 

-4.8180 

2 

.14645 

.44236 

1.6158 

.8652 

-8.928 

.71478 

3.1691 

-.4062 

3 

.30866 

.49990 

1.2906 

-.0235 

-3.636 

.64517 

1.8200 

.0028 

4 

.50000 

.44051 

.9053 

-.5445 

-1.939 

.39879 

.5719 

-.1340 

5 

.69134 

.30843 

.4909 

-.8063 

- 1.048 

.15142 

-.3081 

-.3063 

6 

.85355 

.16199 

.0824 

-1.0163 

-i. 57 it 

.01335 

-.6913 

-.4298 

7 

.96194 

.04499 

-.5951 

-1.0879 

1.643 

-.02677 

-.9888 

-.3230 


The accuracy of this solution with only seven pivotal points Is 
Indicated by comparison with the following values, which were obtained 
analytically for Ag3yu and with I5 plTotal points for Aqa/U; 


X 

u-it 

U " u 

Aqz 

U 

0.14645 

.50000 

.85355 

1.6166 

.9003 

.0725 

- 0.4009 

-.1348 

.4239 


It is seen that the solution using seven pivotal points yields ample 
accuracy. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Dec. 1 , 1954 
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APPENDIX A 
NCfTATION 


A,B 

a 


bn 

C(x) 

CaCx) 


Cp 

cns#<3.ns^®ns^ 1 
^'nsjSna^bQs^ J 

cr 


f(z) 

I 

Ki 

Ka 

M 

N 


Q 

Aqi/^qa 


abscissas of leading and trailing edges, respectively 

factor proportioml to angle of attack measured from ideal 
angle 

coefficient of in series for airfoil ordinate 

camber of airfoil 

camber of fictitious airfoil in second-order solution 

surface-pressure coefficient 

first-order surface-pressure coefficient 

second-order increment in siirface-pressure coefficient 

Influence coefficients for calculating velocity, slope, 
and second derivative of airfoil ordinate 

coefficient in trigonometric polynomial approximation 
to C 

analytic function of complex variable 
imaginary part of f(z) on unit circle 
first-order ccm^ressibility factor, i 

(7 + 1 ) M * + 

second-order compressibility factor, ^ 

free-stream Macb number 

number of subdivisions of chord line in numerical integra- 
tion 

surface speed ratio on parabola in subsonic flow 
(See eq_. (32b),) 

flow speed on surface of airfoil 
first- and second-order increments in q 
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R 

T(x) 

TaCx) 

tr 


U 

u 

V 

X 

Xo 


Y(x) 




y 

z 

a 

P 

7 

’'P 

5 

0 

1 } 

A 

M-p 


real part of f(z) on unit circle 
thickness of airfoil 

thickness of fictitious airfoil in second-order solution 

coefficient in trigonometric polynomial approximation 
to T 

free- stream speed 

velocity perturbation parallel to chord line 
velocity perturbation normal to chord line 
abscissa 

abscissa measirred from edge into airfoil 

distance from round edge measured along initial tangent 
to camber line / 

ordinate of airfoil 

ordinates of upper and lower sixrfaces of airfoil, respec- 
tively 

ordinate 

complex variable 
angle of attack 
/l - 

coefficient in numerical calculation of u 
adiabatic exponent of gas 

coefficient in numerical calculation of Y* 
semivertex angle of sharp edge 

* 

polar angle 

angle of airfoil surface to chord line 
terminal angle of camber line to chord line 
coefficient in numerical calculation of Y" 
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V 


P 

T 


Tl 


tP 




factor that is 1 or i (See eq.. (C12).) 

radius of roiind edge 
airfoil thickness ratio 


3-/3 


T for Joukowski airfoil 


— ^ also (in Appendix C) perturbation velocity potential 



also (in Appendix C) 


perturbation stream function 


( )o 

( )m 
( )x 
( ), 

( )t 

( )c 
"( )" 

n 
c )* 
c )’ 


value at zero Mach nimiber 

value at Mach ntimber M 

first-order approximation 

second-order approximation 

ccanponent associated with thickness 

component associated with camber and angle of attack 

formal series approximation 

uniformly valid approximation 

part not involving angle of attack 

derivative 

value at mth pivotal point counted from trailing edge 
value at nth pivotal point coimted from leading edge 
index of summation, counted from trailing edge 
index of summation, coimted from leading edge 
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APPENDIX B 
AIRFOIL INTEGRALS 

Tbe following are the Cauchy principal value-s for < 1: 
1. (f d| = In 

J x-l 1-x 

-1 


2. 

-1 


3 . 



4. 

_£ ^ ^ ^ ■ 

")-l 



5 . 

-1 -1 

l-(-l) 

n 

6 . 

t 1 d| = 0 

4l Jl _|2 (x-l) 


7- 

f ^ d| = 

'^-1 Jx-l^ (x-l) 


8. 

4 d£ = 

-1 Jl-I^ (x-5) 


9 . 

f dl = 

-1 TLP (x-i) ^ 


10 . 

<r d| = - 3 ix {x^ 

-1 .jCFCx-l) ^ 

-1) 



r ■» 


f t 


10 

H 

• 

to 

0 

H 

VO 

• 

H 

CD 

« 


H 

o\ 

• 

H 

VJl 

• 

• 

H 

(jO 

1 ^ 

• 

• 

H 

H 

H J 

H 

H 

H 

P 

H 

H 

M 

H 

H 
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nl I® In ^ / \ 

£ — if - *7 ■ ^ 


<^1 In i±L . 


2V 2x^ - ^ 


,1 la ^ 


, 1 lija ^ 


pi g m - , r^ ^ 1-6 2 r / \^“^1 V’ 1 

/ "x- T- 5EI1— ^ - n [l-(-D JX n=l=j 


f In — 

^1 (x-6) ^yXiT ^- 41 ^ 


42 - -/l+x 


,1 In 


p- — 1-1 

'Zi n/1-6^ (x- 


d| = - 


n/1-§^ (x-|) •^l-X' 


1 


'£. A-!^ (*-i) 

l. -fl^ (x-l) 

pi S=lu^ 

4 ■ <i| = 

-1 ^yl-x^ (x-6) 


1-X2 


-rt (2 + 3t 


l-x2. 


-Jtx [2 + 3t 


.4 T„ 1+6 

.1 & ^ Xl| 


d6 = - 


^/l-X®y 


3t( ^ + 2x^ + a~ ^ 

V ^ 4 l-x^^ 
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32 . 


/ 


-1 

1 g In 


itl 

izL 


1 (x-l) 


= X 


^yl-|2 (x-i) 


n 

a " ^ 




(-l)''(l)(3)...(v) 


^2^(2v+l)(v.«)^ 


roip 



30 


WACA TN 3390 


APPENDIX C 

THE GERMAIN-WATSON-TEWAITES -WEBER METHOD 


The numerical procedure introduced "by Germain (ref. 21) can he 
adapted to give approximately the thin-airfoil velocities on any profile 
in terms of its ordinates at certain fixed points. In the same way the 
airfoil slope and second derivative can he calculated. Thwaites has 
applied Germain's procedure to thin -airfoil theory (ref. 23)^ and Weher 
has systematized the calculation of the slope and surface velocity for 
uncamhered airfoils (ref. 24) . Here we must treat also cambered airfoils 
and find the second derivative. It is convenient to derive all these 
results from Watson's analysis (ref. 22). 

Let f(z) he regular within the unit circle, and on the unit circle 
have the form 


f(ei0) = R(e) + il(e) (Cl) 

(Our R and I are Watson's itr and e.) Then following Germain, Watson 
approximates to R hy the trigonometric polynomial 


R(0) 


N-i 

as Co + ^ (crcoB r9 + trsln r©) + 


CjtCOS N0 


r=i 


(C2) 


which can he made to coincide with R at the 2N equally spaced pivotal 
points 0 * 0m = mn/N. Thus he derives approximate formulas for I (aside 
from a constant), R' , I', /R, and /l in terms of the values of R at 
the pivotal points times fixed Influence coefficients. 

In thin-airfoil theory the complex perturhatlon potential cp + lilr 
is regular outside the unit circle in the absence of circulation. Inver- 
sion shows that this Involves a change in sign of either the real or imag- 
inary part, since f(e“i0) = R - II. Hence (cp,-i|f) or (\|r,cp) may he iden- 
tified with (R,I) . 

In thin-airfoil theory the tangency condition on the perturbation 
stream function is t = -Y, where Y is the airfoil ordinate. There- 
fore, in order to obtain a solution in terms of the airfoil ordinates we 
identify (Y, -<p) with (R,l) . 
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Streamwise Velocity Increment 


Let X run from 0 at the leading edge of the airfoil to 1 at the 
trailing edge, and 


X = i (1 + cos 0) 
2 


(C3) 


Then the streamwise perturhatlon velocity on the airfoil is given by 


u ^ _ h<p/hd _ 2 ^cp 

U “ “ dx/d0 “ ~ Bin 0 ^ 


Cci^) 


Now according to Watson's eguations (lO), (24), and (2j) , in the absence 
of circulation the values of S(p/B0 at the points 0m are 


2N-1 


'm 


= ^ PpYm+p^ ^p = 


p=0 


2 ’ 

0 , 


N(1 - 


1 

cos 0p ) * 


p = 0 

p = even, not 0 
p =s odd 


(C5) 


Now since PaN-p = Pp^ 


2 N -1 

^ Pp^m+p 

p =0 


N-i 


N-i 


= Pm^o + ^ Pp-mYp + ^ Pp+mY^-p .+ PN-m^N 


(C 6 ) 


p=i 


P =1 


Symmetric airfoils .- For a symmetric airfoil YaN-p = “Yp and 
Yo = Yj]- = 0. Then according to equations (c4), {C^) , and (c6) 


N-i ■ 

- V c Y 
P=i 


cmp = 


sin 0 m 


(Pp-m “ Pp+m) 


(C7) 


This form is convenient for calculating the cmp- It is also easily 
shown, using trigonometric identities, that 
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sin 9 ra,^ 


Cmp = < 0, 


4 sin 0p 

w(cos 0rii - cos 0p)^^ 


p - m = 0 

p - m = even, not 0 ^ 
p - m = odd 


(c8) 


which is Weber's result. 

Antisymmetric airfoils .- For a cambered airfoil of zero thickness 
YaU-p = Yp. The ordinates Yq and Yj^ of the leading and trailing edges 
do not vanish in general, although in the present second-order theory the 
axis is chosen so that they do. Equation (C0 gives 

2N-1 

^ 3pYm+p 

p=0 


= 0m^o 


N-i 

I 

p=i 


(Pp-m + Pp+m)Yp + PN-m^N 


(C9) 


This expression represents the velocity on the \anit circle into which 
the airfoil is mapped. The Kutta condition .will be violated at the trail- 
ing edge of the airfoil imless the expression happens to vanish for m » 0 
Adding a component of circulatory flow changes the velocity on the circle 
by a constant. Hence the Kutta condition is enforced by subtracting from 
the expression of equation (C9) its value at m = 0, so that 

2 H-1 W-1 

Pp^m+p " ^o^^o '*■ 

p=0 


r 

L 

p=i 


^Pp-ra Pp+m 


20p)Yp + On-: 


m - 

(CIO) 


Hence, according to equations (c4) and (C5) 


nm 

U 


N 

= ^ ^mpYp^ 


^mp = 


sin 0m 


p=0 


(Pm ~ Po) f 
(Pp-ra + Pp+m 
(PK-m - Pf)> 


2Pp) > 


p = 0 

p ^ 0, N 
p = N 


? (Cll) 


The expression for dmp can be written more concisely in a form suitable 
for computation as 
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= 


2v 

sin 9m 


(Pp-m + Pp+m ~ ^p)j 


V 



1 otherwise 


n1 


( 012 ) 


Slope of Airfoil 


The airfoil slope is given hy 


Y* = 


dx sin 9 d9 


(013) 


Now according to Watson's equations (29) ^ (3l)» and (3^) 


aN-i 


§1 ■ I 

p=0 


0, p = 0 


(0l4) 


or, since raN-p = “ 7p 


W-i N 


p=0 


P=i 


(015) 


Symmetric airfoils .- Using the symmetry conditions again gives for 
symmetric airfoils 


N-i 


Y^ = ^ ejnpYp, emp = - (7p-m + 7p+m) 


(016) 


p=i 


Antisymmetric airfoils.- Similarly, for camher lines 


N 


Ym = ^ fmpYp, fmp = (7p+m “ 7p-m) 


(017) 


p=0 
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Second Derivative 


The second derivative of the airfoil ordinate is given hy 


V" = ^ 4 ( ^dgy ^ 

^ - dx2 = l^d02 " d0 


(Cl 8 ) 


An approximation for d^/d 02 ^ which is required here, is found by extend- 
ing Watson’s analysis for the first derivative, as he suggests. Following 
closely his section 2.4 gives, after some computation. 



aN-i 

PpYm+p, 

p=0 


tJ-p 


(W - 1)(2N - 1) 

6 


(-DM^n - 


1 - cos 0 p^ 


p = 0 


i 0 


(C19) 


Symm'etric airfoils .- By the foregoing procedure, it is found that 
for symmetric airfoils 

N 

Ym = SmpYp, gmp = [(p-p-m ~ M-p+m) “ cot 0m(7p-m + 7p+m) 1 

p=i (C20) 


Antisymmetric airfoils .- Similarly, for camber lines 
N 

hmpYp, hmp = — [(pp-m + Pp+m) “ cot 0m(7p-m - 7p+m) 1 
P=0 (C21) 


Tables 


The six sets of Influence coefficients required for calculating 
u/U, Y', and Y" for both symmetric and antisymmetric airfoils are tabu- 
lated in tables II and III for IT * 8 . In addition, the coefficients 
required for calculating u/U are given in table V for IT = l 6 . These 
values have been checked by applying them to a number of simple shapes - 
flat plate , ellipse , etc ., - for which the approximation of the airfoil by 
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a trigonometric polynomial is exact. The values are "believed to he accu- 
rate to vithin one unit in the last place. 

For convenience of computation, the coefficients have been renumbered 
so that the pivotal points are counted from the leading to the trailing 
edge. This renumbering is indicated by using indices (n, s) rather than 
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TABLE I.- PIVOTAL POINTS AND ANGLE-OF-ATTACK SOLUTION 



N = 

8 

N = 16 



/1-Xn 



1 

H 

n 

xn 

V Xn 

n 

Xn 

^ — 

0 

0 

00 

0 

0 

CO 

1 

.038060 

5.027339 

1 

.009607 

10.15318 

2 

.i46M^7 

2 . 4 i 4214 

2 

.038060 

5.027339 

3 

.308658 

1.496606 

3 

.084265 

3.296558 

k 

.500000 

1.000000 

4 

.146447 

2 . 4 i 42 i 4 

5 

.691342 

.668179 

5 

.222215 

1.870868 

6 

.853553 

. 4 i 42 i '4 

6 

.308658 

1.496606 

7 

.961940 

.198912 

7 

.402455 

1.218504 

8 

1.000000 

0 

8 

.500000 

1.000000 




9 

.597545 

.820679 




10 

.691342 

.668179 




11 

.777785 

.532511 




12 

.853553 

.414214 




13 

.915735 

.303347 




l 4 

.961940 

.198912 




15 

.990393 

.098491 




16 

1.000000 

0 
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TABHE II.- IKFLUENCE COEFFICIENTS : 


Cns 



1 

2 

3 

4 

1 

20.90501 

-4.07193 

0 

-0.22417 

2 

-7.52395 

11.31371 

-3.35916 

0 

3 

0 

-4.38896 

8 .65914 

-3.15^32 

4 

-.58579 

0 

-3.41421 

8.00000 

5 

0 

- .38896 

0 

-3.15*^32 

6 

-.13291 

0 

-.29769 

0 

n 

0 

-.07193 

0 

-.22417 


®ns 



1 

2 

3 

4 

1 

2 

3 

4 

5 

6 

7 

-6.30865 

17.04789 

-8.92177 

5.65685 

-3.69552 

2.26582 

-1.08239 

-4.99321 

-1.41421 

8.05468 

-4.00000 

2.39782 

-1.41421 

.66364 

1.53073 

-4.71832 

-.44834 

5.65685 

-2.61313 

1.40461 

-.63405 

-0.82843 

2.00000 

-4.82842 

0 

4.82842 

-2.00000 

.82843 

8ns 

i 

1 

2 

3 

4 

D 

-414.3920 

99M83 

46.6274 

-46.8824 

35.3137 
-23.0316 

11.3137 

106.2602 

-160.0000 

76.5287 

-11.3137 

2.0190 

0 

-.2495 

-12.6863 

62.4057 

-97.6081 

54.0559 

-11.3137 

3.8960 

-1.3726 

3.5867 

-11.3137 

50.4692 

-84.0000 

50.4692 

-11.3137 

3.5867 
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FOR THICKNESS, N = 8. 


5 


6 


7 


- .29769 

Alh 21 

.6591^ 

.35916 


-0.07193 

0 

-.38896 

0 

-i)-. 38896 
11.31371 
07193 


. 0 

-.13291 

0 

-.58578 

0 

-7.52395 

20.90501 


5 


0.63405 

-1.40461 

2.61313 

-5.65685 

.44834 

4.71832 

-1.53073 


6 


-0.66364 

1.41421 

-2.39782 

4.00000 

-8.05468 

1.41421 

4.99321 


7 


1.08239 

-2.26582 

3.69552 

-5.65685 

8.92177 

17.04789 

6.30865 


5 


-1.3726 

3.8960 

•11.3137 

54.0559 

■97.6081 

62.4057 

■12.6863 


6 


-0.2495 

0 

2.0190 

-11.3137 

76.5287 

-160.0000 

106.2602 


7 


11.3137 
-23.0316 

35.3137 

-46.882i^ 

46.6274 

99.^783 

-414.3920 
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TABLE III.- INFLUENCE COEFFICIENTS FOR CAMBER, N = 8 



2 

3 

0 

-4.84984 

11.31371 

-4.38896 

0 

.38896 

0 

8.84984 

-11.31371 

-0.43835 

.28130 

-3.75057 

9.05055 

-3.69552 

.87669 

-.57900 

7.10973 

-8.85484 


-.32589 

0 

-3.05260 

8.00000 

-2.60426 

0 

5.98275 

- 8.00000 


5 

6 

-0.19570 

0 

.28130 

-.07193 

-.57900 

0 

.391^1 

-.24509 

-3.69552 

0 

9.53583 

-3.75490 

-3.75057 

11.31371 

7.10973 

4.07193 

-9.097^ ■ 

■11.31371 





3.41421 

-9.22625 

1.41421 

6.16478 

-2.82843 

1.83522 

-1.41421 

1.22625 

-.58579 


-1.61991 

3.69552 

-6.16478 

.4483ii 

5.22625 

-2.61313 



1.00000 

-2.16478 

2.82843 

-5.22625 

0 

5.22625 
-2.82843 
2. 16478 I 
- 1.00000 


-0.72323 
1.53073 

-1.83522 

2.61313 
-5.22625 
-.41)834 
6.16478 
-3.69552 
I.61991I -3. 


415.3214 

-541.5880 

68.1481 

144.5686 

-181.8234 

19^.5097 

-200.0304 

202.5097 

-101.6159 


-4.9706 

80.1516 

- 120.0000 

29.4458 

40.0000 
-52.0732 

56.0000 
-57.5242 

28.9706 


-9.7012 

13.4903 
32.9909 

-66.4121 

21.8234 

21.4903 
-29.1087 

31.^315 

-16.0044 


12.0000 
-22.6274 

16.0000 
22.6274 

-56.0000 

22.6274 

16.0000 
-22.6274 

12.0000 


5 

6 

-16.0044 

31.^^315 

-29.1087 

21.4903 

21.8234 

-66.4121 

32.9909 

13.^903 

-9.7012 

28.9706 

-57.5242 

56.0000 

-52.0732 

40.0000 

29.4458 

-120.0000 

80.1516 

-4.9706 


NACA 
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TABLE rV.- COMPRESSIBILITY FACTORS (7 = T/5) 


M 

Ki 

Ka 

Ka“l 

Ka-l 

2 

0 

1.00000 

1.00000 

0 

0 

.10 

1.00504 

1.01016 

.010162 

.005081 

.20 

1.02062 

1.04271 

.042708 

.021354 

.25 

1.03280 

1.06933 

• .069333 

.034667 

.30 

1.04828 

1.10477 

.104770 

.052385 

.35 

1.06752 

1.15129 

.151294 

.075647 

.40 

1.09109 

1,21224 

.212245 

.106122 

.45 

1.11978 

1.29260 

.292603 

.146302 

.50 

1.15470 

1.40000 

.400000 

.200000 

.55 

1.19737 

1.54654 

.546545 

.273272 

.60 

1.25000 

1.75234 

.752344 

.376172 

.65 

1.31590 

2.05275 

1.05275 

.526373 

.70 

1.40028 

2.51465 

1.51465 

.757324 

.75 

1.51186 

3.27755 

2.27755 

1.13878 

.80 

1.66667 

4.67407 

3.67407 

1.83704 

.85 

1.89832 

7.67085 

6.67085 

3.33543 

.90 

2.29416 

16.1679 

15.1679 

7.58393 

.95 

3.20256 

61.6651 

60.6651 

30.3325 




TAH£ V.- INFLUENCE COEFFICIENTS, N = l6. 


Cub 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


12 

13 

14 

15 

1 

82.(31330 

-15.06135 

0 

-0.63117 

0 

-0.13634 

0 

-0.05070 

0 

-0.02624 

0 

-0.01718 

0 

-0,01344 

0 

2 

-29.5't391 

kl.81002 

-11.20322 

0 

-.70527 

0 

-.18013 

0 

-.07641 

0 

-.04371 

0 

-.03105 

0 

-.02637 

3 

0 

-16.26456 

26.79934 

-8.96045 

0 

-68960 

0 

-.20090 

0 

-.09422 

0 

-.05867 

0 

-.04508 

0 

k 

-2.36017 

0 

-11.42995 

22.62742 

-7.69823 

0 

-.67431 

0 

-.21718 

0 

-.U086 

0 

-.07468 

0 

-.06205 

5 

0 

-1.33236 

0 

-9.05216 

19.24304 

-6.95446 

0 

-.67346 

0 

-.23613 

0 

-.13038 

0 

-.09497 

0 

6 

-.6^366 

0 

-1.14677 

b 

-7-78738 

17.31827 

-6.563P 

0 

-.69189 

0 

-,26237 

0 

-.15668 

0 

-.12424 

7 

0 

-.46165 

0 

-.93529 

0 

-6.96756 

16.31346 

-6.44233 

0 

-.73451 

0 

-.30124 

. 0 

-.19583 

0 

e 

-.85969 

0 

-.36162 

0 

-.80996 

0 

-6.56854 

16.00000 

-6.56854 

0 

-.80996 

0 

-.36162 

0 

-.25909 

9 

0 

-.19583 

0 

-.30124 

0 

-73451 

0 

-6,44233 

16,31346 

^>.967*36 

0 

-.93589 

0 

-.46165 

0 

10 

-.i2Jteii 

0 

-.15666 

0 

-.26237 

0 

-.69189 

0 

-6.56330 

17.31^ 

-7.78738 

0 

-1.14677 

0 

-.64566 

11 

0 

-.09497 

0 

-.13038 

0 

-.33613 

0 

-.67346 

0 

-6.95446 

19.84304 

-9.05216 

0 

-1.53236 

0 

12 

-.06203 

0 

-.07466 

0 

-.11068 

0 

-.21718 

0 

-.67431 

0 

-7.69623 

22.62742 

-11.42995 

0 

-2.36017 

13 

0 

-.04308 

0 

-.05867 

0 

-.09422 

0 

-.20090 

0 

-.68960 

0 

-8.98045 

28.79924 

-16.26456 

0 

Ik 

-.02637 

0 

-.03105 

0 

-.04371 

0 

-.07641 

0 

-.i8(a.3 

0 

-.70527 

0 

-n.20322 

41.81002 

-29.54391 

13 

0 

-.01344 

D 

-.01712 

0 

-.02624 

0 

-.05070 

0 

-.13634 

0 

-.65U7 

0 

-15.06135 

82.01330 


dan 



1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

11 

12 

13 

14 

15 

0 

-66.691^» 

0 

-2.67007 

0 

-0.67654 

0 

-0.31668 

0 

-0.21389 

0 

-0.19329 

0 

-0.24570 

0 

-0.64694 

1 

82.66025 

-18,60789 

.22718 

-1.26820 

*15n!9 

-*33591 

.12868 

-1336B 

.12^ 

-.06350 

-15179 

-.03167 

.22TLB 

-.0^344 

.64694 

2 

-37.14763 

41.81002 

-12.21573 

0 

-1.07882 

0 

-39341 

0 

-.24027 

0 

-.B07M 

0 

-.25823 

0 

-.67332 

3 

.69969 

-17.37779 

£9,044^ 

-9.22666 

.16417 

-.76828 

.13918 

-.22511 

.13916 

-.09422 

.16417 

-.04412 

.84570 

-.01819 

.69969 

4 

-5.84355 

0 

-11.98900 

2.62742 

-7.9498B 

0 

-.83817 

0 

-.35636 

0 

-.26267 

0 

-.30185 

0 

-.76174 

5 

.82379 

-1.92402 

.28928 

-9.1S071 

19*43633 

-6.95446 

.16386 

-.64924 

.163^ 

-.19676 

.19329 

-.08159 

.26926 

-.03170 

.02379 

6 

-2.83771 

0 

-1.52330 

0 

-7.92067 

17.31827 

-6.70247 

0 

-.82058 

0 

-,4l4l6 

0 

-.40238 

0 

-.94803 

7 

1.07227 

-.46165 

.37^3 

-.86673 

*25159 

-6.8^ 

16.58674 

-6.35935 

.21329 

-.64470 

.85159 

-.19644 

•37653 

-.06915 

1.07227 

8 

-1.33216 

0 

-.65089 

0 

-.97413 

0 

-6.65723 

16.00000 

-6.69TO 

0 

-.97413 

0 

-•650ti9 

0 

-1.33216 

9 

1.5^ 

.19583 

•55905 

-.05504 

*37355 

-.53494 

.31668 

-6.25795 

16.63014 

-6.77912 

.37355 

-.70333 

.55905 

-.15666 

1.59205 

10 

-.94803 

0 

-.40238 

0 

-.41416 

0 

-.82058 

0 

-6.70247 

17.31827 

-7.92067 

0 

-1.52330 

0 

-2.23771 

11 

2.88338 

1.01826 

1.01251 

.48665 

.67654 

.23613 

•57»4 

-.24744 

.57354 

-6.51954 

19.91958 

-0.55247 

1.01251 

-.87405 

2.88338 

12 

-.76174 

0 

-.30185 

0 

-.26267 

0 

-.35636 

0 

-.83817 

0 

-7.94982 

22.6^J|2 

-U.989W 

0 

.5.24355 

13 

7.60372 
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r^ure 2.- 'Speed on t2- percent -thick ^mmetricat Joukowskt airfoil at zero 
angle of attack in incompressible flow. 
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Figure 4r Speed on NACA 0012 airfoil at zero angle of attack in incompressible flow. 
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Figure 5.- Pressure coefficient on NACA 0012 airfoil at zero angle 


of attack in subsonic flow. 
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Figure 6.- Comparison of theoretical and experimental pressure distributions 
on NACA 0015 airfoil at M=0.70, zero angle of attack. 
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